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Within an effective nonlinear chiral Lagrangian framework the substructure of /o(1370) is studied. 
The investigation is conducted in the context of a global picture of scalar mesons in which the 
importance of the underlying connections among scalar mesons below and above 1 GeV is recognized 
and implemented. These connections are due to the mixings among various quark-antiquarks, four- 
quarks and glue components and play a central role in understanding the properties of scalar mesons. 
Iterative Monte Carlo simulations are first performed on the 14-dimensional parameter space of the 
model and sets of points in this parameter space (the global sets) that give an overall agreement 
with all experimental data on mass spectrum, various decay widths and decay ratios of all isosinglet 
scalar states below 2 GeV are determined. Then within each global set, subsets that give closest 
agreement for the properties of /o(1370) are studied. Unlike the properties of other isosinglet states 
that show a range of variation within each global set, it is found that there is a clear signal for 
/o(1370) to be predominantly a quark-antiquark state with a substantial ss component, together 
with small remnants of four-quark and glue components. 

PACS numbers: 14.80.Bn, ll.30.Rd, 12.39.Fe 


I. INTRODUCTION 


The internal structure of light scalar mesons has continued to challenge the quark model and QCD for several decades 
[H . The light and inverted mass spectrum for the lowest scalars does not obey a simple quark-antiquark picture which 
is known to work reasonably well for other spins such as pseudoscalars and vectors. The MIT bag model of Jaffe 0, 
which is based on a diquark-antidiquark picture, explains why these states are lighter than expected and provides 
a natural template for their inverted mass spectrum. While the states above 1 GeV are expected to be closer to 
quark-antiquark states, nevertheless they too show deviations from such a simple scenario. Particularly, the case of 
isosinglet scalars are more involved than other isospin channels for they mix with glue in addition to the two- and 
four-quark components. A wide range of investigations on scalar meso ns (t hat include lattic e QCD, Q CD sum-rules, 
chiral models and effective theories) can be found in the literature [^- |l2l| (some of which (l06j| - [l2ll| have /o(1370) 
as one of their focus points). 

Important information on scalars are obtained in various pseudoscalar scatterings such as in tttt channel for studies 
of the isosinglet states [particularly /o(500), or sigma (cr)]; in ttAT channel for the studies of the isodoublet states 
Arg(800) [or kappa (k)] and Aro(1430); and in ttt ] channel for the studies of the isovector states ao(980) and ao(1450). 
Chiral Lagrangians provide an effective framework for investigating pseudoscalar interactions. Particularly, chiral 
perturbation theory provides a systematic approach to studies of pion interactions near threshold. In this ap¬ 
proach, pions are the main fields of interest and therefore the heavier fields of vectors and scalars are integrated out. 
However, for the purpose of exploring the properties of scalar mesons, which are outside the immediate focus of chiral 
perturbation theory, it is natural to explicitly keep the scalar meson fields in the Lagrangian instead of inte grating 
them out. Two suitable frameworks, that are the foundation of the present s tudy , are the linear sigma model |74l - l^ 
as well as the nonlinear chiral Lagrangian models that include scalar fields 9ll - ll01 |. In such model buildings, the main 
guiding principles are the well-known chiral symmetry and its breakdown, isospin symmetry (and in relevant processes 
its breakdown), the U(1)a axial anomaly, and the main assumptions that need to be made are related to modeling the 
QCD vacuum as well as the potential. The choice between linear versus nonlinear is a matter of the processes to be 
investigated and the information to be extracted, nevertheless, they are overall complementary. Prior works by one 
of the authors within the linear and the nonlinear models (some of which listed in refs. have indeed shown 

a consistent pattern for the scalar mesons. Specifically, the properties of sigma and kappa extracted in the nonlinear 
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model in m and [ 9 ^ respectively, are quantitatively close to those found within the linear sig ma model [s^; the 


aatterns among the quark-antiquark and the 
I are consistent with similar patterns studied 


four-quark nature of light scalar mesons below 1 GeV studied within the nonlinear model in |93l | are consistent with 
the results within the linear model in [s^; and the underlying mixing 
four-quark components of scalars below and above 1 GeV studied in [98 
in detail within the generalized linear sigma model of refs. [82l - [8^ . 

The approach we take in this paper is within the nonlinear model. In the context of the nonlinear chiral Lagrangian 
of refs. several low-energy processes that probe light scalar mesons are investigated. In order to describe 

the experimental data within this framework, there is a need for a cr and a k in the analysis of tttt M and ttK 
scattering, respectively. Motivated by the evidence for a cr and a ft, and taking into account other experimentally 
well-established scalars [the /o(980) and the ao(980)] a possible classification of these scalars (all below 1 GeV) into a 
nonet is studied in and it is shown that there exists a unique choice of the free parameters of this model which in 
addition to describing the tttt and ttK scattering amplitudes, well describes the experimental measurements for several 
decays, such as, for example, the rj' —^ rjinr decay [9^. The insight into the quark substructure is obtained through 
the mixing patterns between the properties of the isosinglets [a and /o(989)]. The best value for the isosinglet scalar 
mixing angle found in [ 9 ^ is clearly consistent with an ideally mixed qqqq assignment of the MIT bag model i)- 
Various predictions within the framework of ref. 93| are in close agreement with other experimental or theoretical 
works. These include the 70 MeV estimate of the total decay width of ao(980) in 19411 which is confirmed experimentally 


[65j ; the four-quark nature of light scalars probed in radiative decays 


[; and the prediction of sigma meson 


in agreement with experimental analysis |66l |67| . Other low-energy processes in which scalar mesons are expected 
to play important roles are studied within this framework including the ttt; scattering and the isospin violating 
decays r], q' —>■ Stt [o^. Overall, the framework of (9^ has resulted in a coherent description for the physics of scalar 
mesons below 1 GeV. 

The first step in extending the framework of to include the scalars above 1 GeV (in addition to those below 1 
GeV) was done in ref. in which two scalar meson nonets (a quark-antiquark nonet and a four-quark nonet) were 
introduced and the properties of / = 1/2 and I — 1 scalar mesons based on an underlying mixing among the two- and 
four-quark nonets were studied. The idea that was introduced in is rather simple, but quite effective: Assuming 
that there is a four-qaurk scalar meson nonet below 1 GeV (which was proposed in MIT bag model 0 and has been 
since supported by many independent works), as well as a quark-antiquark scalar meson nonet above 1 GeV (which is 
expected to be a reasonable template for some of the scalars above 1 GeV), then it is natural to investigate whether 
some of the observed deviations in mass and decay widths of the scalars above 1 GeV can be related to a mixing 
among these two- and four-quark nonets. It was shown in how such a mixing leads to a “level-repulsion” that 
explains: (a) why scalars below 1 GeV are so light, and (b) why there are unexpected deviations in mass spectrum 
and decay widths of the I = 1/2,1 scalars above 1 GeV (a brief review of this mixing mechanism is given in the 
Appendix). 


Further steps in extending the framework of [98| to include the 1 = 0 states above 1 GeV were taken in 
in which various additional mixi ngs with glue are also present, and as a result, make the analysis considerably more 


9i-[InI| 


challenging. Specifically, in |99l . |100| a preliminary study of the mass spectrum as well as various decay widths and 
decay ratios of isosinglet scalars were given, however, for simplicity, an uncorrelated analysis was performed (the 
mass matrix for this system is a 5 x 5 matrix of eight a priori unknown parameters and the scalar-pseudoscalar- 
pseudoscalar vertices depend on six coupling parameters, therefore at first it may seem that the mass spectrum and 
the decay analyses are uncorrelated, but since the physical states that appear in the decay analysis are obtained by 
diagonalizing the mass matrix, the decay analysis implicitly depends on the mass matr ix parameters and consequently 
this establishes a correlation between these two sets of parameters). In the work of |l0l| the effect of large experimental 
uncertainties on some of the scalar masses on determination of quark and glue components were studied in detail. 
In the present work, we use an iterative Monte Carlo method (developed by the authors) to extend the works of 


99l llOPjl to a correlated analysis of the mass spectrum, decay widths and decay ratios of all / = 0 scalars below 2 GeV 
and simultaneously examine the 14 unknown parameters in the mass and interaction parts of the Lagrangian. This 
results in estimating the substructure of all / = 0 scalar states below 2 GeV, which are, in general, rather sensitive 
to the experimental inputs and the physical conditions imposed. However, among the five isosinglet states [/o(500), 
/o(980), /o(1370), /o(1500) and /o(1710)], the exception is /o(1370) which is found in our simulations to exhibit the 
mos t sta ble substructure and to be dominantly a quark-antiquark (mainly ss) state, further supporting prior findings 

After a brief review of the model in Sec. II, we set up the numerical work in Sec. Ill followed by results in Sec. IV 
and a summary in Sec. V. 
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II. BRIEF REVIEW OF THE MODEL 


The scalar properties are typically probed in various pseudoscalar scatterings or in their decay channels to pseu¬ 
doscalars. Therefore the light pseudoscalars (tt, K and rj/t]') are essential ingredients in models for investigating scalar 
mesons, and the model we are using in this work is no exception. It will contain the psedudoscalars below 1 GeV 
as well as scalar mesons below and above 1 GeV. In certain processes such as tttt scattering, the vector mesons also 
contribute and in theose cases vectors are added too. The leading pseudoscalar Lagrangian density is (see [^ 1 

£0 = -^Tr {d^Ud^U^) + Tr [B {U + U^)] , (1) 


with 


U = 


( 2 ) 


where (j) is the conventional matrix of pseudoscalar fields and = 0.131 GeV is the pion decay constant. The second 
term is the symmetry breaking term with B = diag{Bi, Bi, B^), with Bi = and B^ = F^{mj^ — m^/2)/4. 

Moreover, the U(1)a breaking terms induced by instanton effects need to be added to the Lagrangian density ([T|) to 
generate the rj' mass 


£ 


v' ~ 


c 


In 


/ detU 
\detU^) 


(3) 


where c is a constant proportional to the rj' mass (the dots represent additional terms given in Eq. (2.12) of [H). 
Note that the functional form of Eq. expressed in terms of In and det functions, schematically shows that chiral 
SU(3 )lX SU(3)r symmetry is maintained while the U(1)a is broken to generate a mass term for ry'. 

Under chiral transformation 


U UlUUI- 


(4) 


The nonlinear model is derived by integrating out the heavy fields of scalar mesons and gives a convenient framework 
for investigating the Goldstone bosons interactions near the threshold. However, for our present objective of exploring 
the properties of scalar mesons we need to reintroduce the scalars back into the Lagrangian. One way of course is 
to start from the linear sigma model, form chiral nonets and study the interactions. This is done in refs. [ziQ 
and involves the heavy pseudoscalars 7r(1300), iL(1460) and several etas [77(1290), 77(1405) and ?7(1760)] which are 
related to the chiral partners of heavy scalar mesons. An advantage of the nonlinear approach is that it simplifies the 
framework by not bringing in the heavy pseudoscalars which are actually not needed in the mass spectrum and decay 
analyses of the present work. The disadvantage of the nonlinear model is that it looses contact with the “nuts and 
bolts” of chiral symmetry and its spontaneous breakdown via QCD vacuum condensates, which is so transparently 
traced in the linear model. Nevertheless, for the purpose of analyzing the mass spectrum and decay channels of 
isosinglet states, the no nlin ear model is reasonably effective. As stated before, in general, various analyses within the 
nonlinear model in [9ll - ll01 | have been consistent and complementary to those in the linear model [7414^ . 

It is noted that under chiral transformation the field ^ defined by C/ = transforms as 


e ^ Ul^IC^ (<)., Ul, Un) = 1C (<^, Ul, Ur) (Ul 


(5) 


where K, is defined in the above equation. Then it is easy to show that under chiral transformation the object 


£ 7 . = ^ , (6) 

transforms as 

Pfj, ^ ICpf^K.!. (7) 

To reintroduce the scalars back into the nonlinear Lagrangian, it was considered in that the scalar nonets were 
made of “constituent” quarks and transform in the same way as © 

n^k:nk:\ 


( 8 ) 
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where TV is a four-quark scalar meson nonet [93l| which is defined in terms of diquark-antidiquark fields 


according to 


Qa — ^abcQ Q ; 



abc 

— 6 

qtqc, 



dsds 

dssu 

dsud 

oc QaQ'" = 

suds 

susu 

suud 


udds 

udsu 

udud 


(9) 


( 10 ) 


The transformations 0 and (|8]) allow writing the Lagrangian terms for scalar fields as well as the scalar-pseudoscalar- 
pseudoscalar interaction terms which involve nonlinear pion terms (for details see Appendix B of (9^). Similarly, a 
conventional quark-antiquark nonet N'^ oc (that has exactly the same transformation property, and hence can 
mix with TV) can be added 


N'^ICN'K). (11) 

This was introduced in ref. (9^ in which the case of / = 1/2,1 states were studied in detail. 

Here, the relevant terms for our analysis are the scalar mass terms and the scalar-pseudoscalar-pseudoscalar inter¬ 
action terms which are extracted from the general chiral invariant Lagrangian of ref. together with the mixing 
mechanism of [9^ . In addition a scalar field G that represents the effective field of a scalar glueball and is relevant to 
th e stu dy of isosinglet states is also added [9^. Identification of G with scalar glueball, which is discussed in detail 
in |l24| . is based on the results of various fits to experiment which show that, for example: (i) /o(1500) and /o(1710) 
are the only two states that contain a high content of this field; (ii) G couples strongly to pseudoscalar-pseudoscalar 
channels that involve rj'] and (iii) the mass of G is determined to be in the range of 1.5-1.8 GeV consistent with lattice 
QCD estimates. These observations are all consistent with identifying field G with an scalar glueball. The mass and 
interaction parts of the Lagrangian are 


C 


1=0 

mass 


/’i=U _ 


_ cTr(TV)Tr(TV) - dTr(TV)Tr(TVAl) 

-c'Tr(TV')Tr(A^') - d'Tr{N')Tr{N'M) - gG^ - pTr(A)Tr(A^') - eGTr (TV) - /GTr (TV'), (12) 

+ BTt (TV) Tr (9^(/9^<).) + DTy (TV) Tr (5^</.) Tr (9^</.) + H'Tr (TV') Tr 
+D'Tv (TV') Tr (5^</) Tr (5;,</) + AGTr (5^</5;,</) + FGTr (5^,/) Tr (5^,/), (13) 


,/=l/2, 

.'mass 


where the matrix Ai = diag(l, l,x) with x being the ratio of the strange to non-strange quark masses, and Cn 
and are the mass terms and the scalar-pseudoscalar-pseudoscalar interaction terms relevant to the I = 1/2,1 

states studied in ref. [9^ that also contribute to the 1 = 0 states. Therefore, part of the Lagrangian of isosinglet states 
is constrained by the properties of / = 1/2,1 scalar mesons. All free parameters in J = 1/2,1 parts are determined 
in fits to the mass spectrum and decay properties of / = 1/2,1 states in (these parts are briefly presented in 
Appendix A). It is also shown in the Appendix A that the mixing term between the two- and the four-quark nonets 
TV and TV' is similar to the instanton contribution to the scalar sector which is studied in the literature and 

is suggested to be important for the isosinglet scalar states. We see that the investigation of isosinglet states is not 
independent of the isodoublets and isovectors and indeed is constrained by them. The rest of the Lagrangian density 
only contribute to the 1 = 0 states and is considerably more complex than the I = 1/2,1 parts due to the fact that 
there are internal mixings of two isospin zero combinations within each nonet TV and TV', as well as the mixings of 
these combinations with a scalar glueball, and as a result, there are more parameters to keep track of and the mixing 
matarix is 5 x 5 as opposed to 2 x 2 mixings in the case of isosinglets and isovectors. Initial studies of the 1 = 0 
sector are given in 9^101 | and will be generalized in the present investigation. There are altogether 14 unknown 
parameters, eight of these (c, d, c', d', g, p, e and /) only contribute to the isosinglet 5x5 mass matirx (in addition, 
there are also contributions to the 1 = 0 mass matrix coming from jCmass), and the remaining six parameters (B, D, 
B', D', E and F) contribute to the scalar-pseudoscalar-pseudoscalar coupling constants (there are also contributions 
to the interaction vertices coming from 

Our convention for the mass matrix is as follows: Altogether, there are five 1 = 0 combinations, two in the four- 
quark nonet [TVg and (TV/ -|- TV|)/\/2], two in the quark-antiquark nonet [TV'g and (TV'/ -|- TV'2)/\/2] and a scalar 
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glueball (G). We organize these components into a column matrix 


Fo = 


/ Nl \ 

{Nl + NI)/V2 
N’l 

{N'l+N'l)/V2 

\ G J 


/ \ 


/ 

udud 

\ 

/o" 



{dsds + susu) j \p2. 


fis 

Jo 

(X 


ss 


f/NS 

Jo 



{uu + dd) / \/2 


\ G J 


1 

asG^^^G^, 

/ 


(14) 


where the corresponding quark substructures are shown on the right and S and NS stand for strange and non-strange. 
Then the mass terms are organized into the mass matrix 


- ifoM^Fo = -ipMLgF, 


where F contains the five isosinglet physical fields 


F = 


/ /o(500) \ 
/o(980) 
/o(1370) 
/o(1500) 

V /o(mo) / 


= K-^Fo, 


(15) 


(16) 


and Kp^ is the rotation matrix that converts the quark and glue basis into the physical basis. 

The unknown parameters c and d induce “internal” mixing between the two 1 = 0 flavor combinations [{Nf+N^)/ 
and N^] of nonet N. Similarly, c' and d' play the same role in nonet N'. Parameters c, d, c' and d' do not contribute 
to the mass spectrum of the I = ll2 and 7=1 states. The “external” mixing between nonets N and N' (the p term), 
the glueball mass term (the g term), and the glueball mixing terms with nonets N and N' (the e and / terms) are 
also given in Eq. m- Parameters B and D are unknown coupling constants describing the coupling of the four-quark 
nonet N to the pseudoscalars, parameters B' and D' are couplings of N' to the pseudoscalars, and parameters E and 
F describe the coupling of a scalar glueball to the pseudoscalar mesons. 

After diagonalization of the mass matrix and rotation of the quark and glue basis to the physical basis, the 
interaction Lagrangian (nsi) can be rewritten as: 

-Ant = + 

+%n F,dpridf,ri + Fid^gd^g' + F.d^g'd^g', (17) 

where y*g, is the coupling of the i-th isosinglet scalar [with i = 1... 5 corresponding to the physical states 
/o(500)... /o(1710); see Eq. (ITBl) ] to pseudoscalars s and s', and is given by 


j 

with K defined in (HU) and -fss' = diag (y,^,^, yf,,, yf,, , y£,) , in which the diagonal elements are the couplings 

of the pseudoscalars s and s' to the quark and glue basis , fo , f'o, f'o^: G [defined in (fT4l) ]. respectively. The 
diagonal elements for all decay channels ss' are given in . In next two sections we give the details of our numerical 
determination of the 14 free parameters in the Lagrangian density. 


III. SETTING UP THE NUMERICAL ANALYSIS 


There are 14 unknown parameters in the 7 = 0 part of the Lagrangian density that we need to determine by 
incorporating appropriate experimental data on the mass spectrum as well as the appropriate decay widths and decay 
ratios. These can be divided into a six-dimensional parameter space (77, D, B', D', E and F) that only affect the 
scalar-pseudoscalar-pseudoscalar coupling constants, and an eight-dimensional parameter space (c, d, c', d', g, p, e 
and /) that both directly enter into the 5x5 mass matrix, as well as indirectly enter in the calculation of decay 
widths and decay ratios through the rotation matrix K that rotates the “bare” bases into the physical bases. As a 
result, determining these two groups of parameters independent of each other is on ly an approximation and the exact 
determination requires a simultaneous 14-parameter fit. In works of refs. HiMl, as a preliminary approach, these 
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Short notation 

Quantity 

Experimental value [ref.] 

mi 

m[/o(500)] 

400 - 550 MeV [IJ 

m2 

m[/o(980)l 

990 ± 20 MeV [1] 

m 3 

m[/o(1370)] 

1312 MeV [122] 

7724 

m[/o(1500)] 

1502 MeV [122] 

7725 

m[/o(1710)] 

1727 MeV [12^ 

KK 

r[/o(1370)-»7r7rl 

rIfn(1370)-»KKl 

2.17 ±0.9 [12^ 

p3 
i vv 

KK 

r[/o(1370)^7,7,l 

r[/o(1370)-»KKl 

0.35 ± 0.30 [12^ 

rA, 

nv 

rf/o(1500)-^7r7rl 

r[/o(1500)^i77)] 

5.56 ±0.93 [122] 

p4 

J- KK 

r[/o(1500)-)-K/C] 

r[/o(1500)-^7r7r] 

0.33 ± 0.07 [122] 

r4 

^ njii. 

r[/o(1500)^i77)'] 

rl/o(1500)^i77)l 

0.53 ±0.23 [122] 

KK 

r[/o(1710)-»vr7r] 

rIfn(1710)-»KKl 

0.20 ± 0.03 [122] 

p5 

r[/o(1710)->.7i7l 

r[/o(1710)^/fKl 

0.48 ±0.19 [12^ 

rL 

r[/o(500) — >■ tttt] 

400 - 700 MeV [1] 

rL 

r[/o(980) ->• tttt] 

40 - 100 MeV yj 

r3 

TTTT 

r[/o(1370) ^ tttt] 

(0.26 ± 0.09) X (230 ± 15) MeV [116] 

KK 

r[/o(1370) ^ KK] 

(0.35 ± 0.13) X (230 ± 15) MeV [116] 

r4 

TTTT 

r[/o(1500) —>■ tttt] 

(0.349 ± 0.023) X (109 ± 7) MeV [1] 

p4 

J- KK 

r[/o(1500) KK] 

(0.086 ± 0.010) X (109 ± 7) MeV [1] 

r4 

riri 

r[/o(1500) rj-ij] 

(0.051 ± 0.009) X (109 ± 7) MeV [1] 

r4 

^ 7777 / 

r[/o(i50o) ^ w/] 

(0.019 ± 0.008) X (109 ± 7)MeV[lJ 

r® 

TTTT 

r[/o(1710) ->• tttt] 

(0.12 ± 0.11) X (220 ± 40) MeV [123] 

pS 

J- KK 

r[/o(1710) ^ KK] 

(0.36 ± 0.12) X (220 ± 40) MeV [123] 

r® 

r[/o(i7io) ^ w] 

(0.22 ± 0.12) X (220 ± 40) MeV [12^ 


TABLE I: Target quantities used to explore the 14 parameters of the Lagrangian in global fits I and II. In global fit I (II) the 
decay channels of /o(1370) are exclnded (included). The short notation for the quantities are defined in column one. 


two parameter spaces were studied in separate fits in some details. Here we generalize those separate fits into one 
simultaneous fit using an iterative Monte Carlo algorithm. 

Since the experimental status of scalar mesons is not yet firmly established, different existing data do not always 
overlap. To reduce uncertainties that stem from unestablished experimental data, we incorporate several sources of 
input. We have collected the inputs that we use into two Tables U and m The experimental inputs in Table U are 
divided into three gronps: the masses and several decay ratios and decay widths. Altogether, 23 inputs are displayed 
in Table HI In principle one might think about selecting a subset 14 out of these 23 experimental data and examine 
whether a determined system of 14 equations in 14 unknowns might be formed. However, we do not seek to solve a 
mathematical system, because even if such a determined system exists, solving such a highly nonlinear system and 
finding all distinct solutions, can be at the expense of pushing the model predictions away from other experimental 
quantities not included in the chosen set of 14 inputs. Our objective is to explore the underlying mixings among 
various two-quark, four-quark and glue components, in order to achieve a global understanding of all / = 0 states. 
This objective is sometimes at the expense of individual accuracies, at least at the present approximation of the 
model. Therefore, we aim to determine the 14 Lagrangian parameters such that we get a reasonable agreement with 
the experimental data displayed in Table HI We use the inputs of Table HI in two ways: In our global fit I, we exclude 
the two partial decay widths of /o(1370) and input the rest of 21 quantities. This is because of the fact that even 
though these two decay widths are listed in PDG [I| they are not used in any averaging and therefore we would like 
to test their importance. In global fit H we include all the 23 quantities (including the two decay widths ignored in 
the global fit I). Our global fit HI is obta ined from the inputs of Table IH] which does not include the decay ratios 
provided by the WA102 collaboration |l22l | displayed in Table HI 
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Short notation 

Quantity 

Experimental value [ref.] 

mi 

m[/o(500)] 

400 - 550 MeV [IJ 

m2 

m[/o(980)l 

990 ± 20 MeV [1] 

m3 

m[/o(1370)] 

1300 ± 15 MeV [116] 

7724 

m[/o(1500)] 

1505 ± 6 MeV [Ij 

7725 

m[/o(1710)] 

1690 ± 20 MeV [12^ 

Pin 

r[/o(500) tttt] 

400 - 700 MeV [Lj 

Pin 

r[/o(980) tttt] 

40 - 100 MeV [IJ 

J .3 

r[/o(1370) tttt] 

(0.26 ± 0.09) X (230 ± 15) MeV [1161 

pS 

KK 

r[/o(1370) ^ KK] 

(0.35 ± 0.13) X (230 ± 15) MeV [116] 

r4 

r[/o(1500) — >■ tttt] 

(0.349 ± 0.023) X (109 ± 7) MeV [Lj 

r4 

^ KK 

r[/o(1500) KK] 

(0.086 ± 0.010) X (109 ± 7) MeV [1] 

r4 
^ 777 ? 

r[/o(1500) gg ] 

(0.051 ± 0.009) X (109 ± 7) MeV [1] 

r4 

^ 7777 / 

r[/o(1500) -s- ggi ] 

(0.019 ± 0.008) X (109 ± 7)MeV[lJ 

r® 

TTTT 

r[/o(1710) ^ tttt] 

(0.12 ± 0.11) X (220 ± 40) MeV [123] 

pS 

KK 

r[/o(1710) ^ KK] 

(0.36 ± 0.12) X (220 ± 40) MeV [123] 

r® 

^ 2777 

r[/o(1710) ^ gg ] 

(0.22 ± 0.12) X (220 ± 40) MeV [1^ 


TABLE II: Target quantities used in global fit III. 


We measure the goodness of the fits by the smallness of the parameter y defined by 


/V“P 

q 

X(pi.. .pii) = ^ 


r" - ipi • ■ - pm) 

.exp 


(19) 


where ± with i = 1... iV^P (for quantities that an experimental range is reported, g™’’ = 

Central value). Our target quantities are q^^^ which are also theoretically 
calculated by the model as a function of the 14 model parameters {pi.. .ps = c, d, c', d', g, p, e, f and pg .. .pi 4 
= B, D, B', D', E and F). This results in a “fixed target method” in which the computation revolves around repro¬ 
ducing the central values of the experimental data, which here in this work is just one set. This is to be contrasted 
with a “moving target method” in which every point within the experimental range is treated as a viable target and 
the computation spans over all possible ta rget sets to find the best agreement with the model computation. The 
guiding function %, that was introduced in |lOl| , has two important advantages that are suitable for our study of the 
underlying mixings and the global picture of scalar mesons: (a) it gives each individual data an equal weight (unlike 
the conventional method); and (b) it can be used when dealing with different quantities, such as here that we 
input different types of experimental quantities of masses, decay widths and decay ratios. To measure the goodness 
of our fits, we compare a x with its corresponding experimental value defined by 




i^exp 

q 

E 

2 = 1 


Ag- 




( 20 ) 


Since in most situations there is nothing unique about the experimental central values, we do not limit our simulations 
to just searches for the best fit (or lowest y value). Instead, we adopt an inclusive process that accounts for the 
experimental uncertainties (expressed by x®’'p) by filtering out simulations that do not satisfy y < y®’^P condition. 
This method in general leads to finding an acceptable set instead of just the best point which obviously is also included 
in the set (the method results in “dispersive fits” that reflect the experimental uncertainties, but for simplicity we 
refer to them as “fits” throughout this work). In the moving target method, for each choice of an experimental set, 
there is a corresponding y®’'^, and therefore, overall this method leads to a range for 

For the details of our numerical analysis we will use several strategies and define the guiding function y accordingly. 
In general our guiding function contains three parts 


y (pi • ■ -pii) = Xm [pi ■ ■ -Ps) + Xr {pi ■ --Pia) +X(r/r) (pi ■ ■■Pia) , 


(21) 




























where the three terms on the right refer to x for mass, decay width and decay ratio, defined by 
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(22) 


(23) 


where i = 1... 5 correspond to the five isosinglet scalars in ascending order of masses, and a and /3 are the two-body 
decay channels and in this work take values 1... 4 which respectively correspond to the decay channels tttt, ■kK, rjr] 
and rjr]' (note: it is understood that the summations run over relevant values of a and /3 that are listed in Tables U 

andHTl). 


IV. RESULTS 

The investigation of the substructure of /o(1370) in the present work, is based on first establishing global rela¬ 
tionships among all relevant scalars below and above 1 GeV and then zooming in on /o(1370). With that in mind, 
we first perform global fits of the model predictions to experimental data for the quantities given in Tables U and HIl 
and search within the 14-dimensional parameter space for sets of acceptable points that satisfy an overall agreement 
between model predictions and experiment. Then within these acceptable sets, we further zoom in on properties of 
/o(1370). This leads to our global fits and their further refinements in this section. 


A. Global fits 

To obtain the acceptable points within our 14d parameter space that provide an overall acceptable description of 
all relevant scalars, we perform three global fits. In global fit I, we exclude the two decay widths of /o(1370) from the 
23 target list of Table |T] and in global fit II we include all the 23 inputs of that Table. Global fit III is obtained with 
the target values of Table [TTl 


1. Global fit I 

In global fit I, where the partial decay widths of /o(1370) are excluded from the target inputs, the guiding function 
for XI is computed from m in which Xm, Xr and X(r/r) are obtained from (12211 . with the condition that in xr the 
decay widths of /o(I370) have been excluded (i.e. i ^ 3). The Xm and X(r/r) are those given in (1^ and include 
all the data given in Table HI We use Monte Garlo simulation over the 14d parameter space and search for points 
p = {pi .. - Pli) for which 

Xiip) < XT^, (24) 

subject to the constraint 

r[/o(1370) ^ (tttt -k KK + r/r/)] < 500 MeV. (25) 

In this case x™^ = 7.3. This leads to a set of points 

^1 = {p\p G : conditions (l24ll and (1251) are upheld} . 


(26) 
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8. Global fit II 

In global fit II, that we include the tttt and KK partial decay widths of /o(1370), the guiding function (yn) is 
computed from (ED with all data in Table U included in Xm: Xr and X(r/r)- For tFis case too, we use the same 
iterative Monte Carlo method to search through the 14-dimensional parameter space for points at which 

Xnip)<Xu^- (27) 

In this case = 8.2. This leads to the second set 

'S'li = {p Ip G : condition (l?7l) is upheld} . (28) 


3. Global fit III 


In global fit III, target inputs are given in Table Ull in which the decay ratios given by WAI02 collaboration 1221 
are not included. Similar to the previous two global fits, we use the same iterative Monte Carlo method and scan the 
14d parameter space for points at which 


XIII (P) < XnT- 


(29) 


where Xin’ = 5.6. This leads to the third set 

‘S'lii = {pIp G : condition (1^ is upheld} . (30) 

For the three global fits I, II and III, the results are displayed in Fig. [T] and compared with their corresponding 
experimental data. Our Monte Carlo simulations show that the points in the 14d parameter space (squares) that 
satisfy the global conditions (l24t and (1251) for fit I, condition ED for fit II, and condition (l29l) for fit III, lead 
to properties of /o(1370) that generally overlap with the experimental data (solid circles and error bars) with the 
exception of the input mass of /o(1370) that in our simulation comes out larger than its target experimental values 
displayed in Tables U and |lll We interpret this deviation of mass from its target values as a measure of the size of 
the next order corrections beyond the present leading order of the model which falls in the range of 8% — 21%. Also 
shown are the simulation averages and one standard deviation around the averages (triangles and error bars). 

Fig. [5] shows the quark and glue components of /o(1370) in the basis defined in (jlT)) . In this figure the individual 
points (shown by “-I-”) are just the results of the global fits without imposing any additional conditions. Also shown 
are the averages (triangles) and standard deviations (error bars). The components are respectively proportional to 
udud, {dsds + suus)/\/2, ss, [uu + dd)/'/2 and glue G. For convenience, the total percentages of the four-quark, the 
quark-antiquark and the glue for each case are also given beneath their detailed component figure in Fig. [2] Clearly, 
the global fits show that the four-quark and glue components are considerably smaller than the quark-antiquark 
components. The ss is the dominant component in either case. We will of course further zoom in on the properties of 
/o(1370) in the next subsection, however, we will see that the conclusion will remain unchanged. For the convenience 
of the reader, we have also made histograms in Fig. |3]that show the distribution of the simulations for each component 
by breaking them down into five bins. This makes it easier to see that the two quark components (particularly ss) 
are the only ones that have high-percentage bins filled. The substructure of /o(I370) shown for the global fits I, II 
and III in Figs. [2] and |3] has been a consistent and stable trend in all our simulations, regardless of the additional 
conditions or filters imposed. This is somewhat in contrast with the case of other isosinglets studied in this model 
|/n^ 00), /o(980), /o(1500) and /o(1710)] for which their substructure are rather sensitive and closely correlated (see 


B. Zooming in on /o(1370) 

In previous subsection we made three global fits to the mass spectrum and decay properties of isosinglet scalar 
states below and above I GeV. This led to identification of three sets of points in the I4d parameter space [S':, S':: and 
S'!!! defined in (E51) . (E51) and (1501) . respectively] for which there is an overall agreement between the model predictions 
and their corresponding experimental data. Working within these three sets ensures that the model respects (at least 
qualitatively) the global properties of all of these isosinglet scalar states. In this subsection, we further refine our 
investigation of /o(1370) by searching for subsets within each of these three sets that better describe this state. 
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FIG. 1: Mass and decay properties (GeV) of /o(1370) obtained from Monte Carlo simnlation for global fits I (left), II (middle) 
and III (right) defined in subsection IV-A. The results of simulations (squares) and their averages and one standard deviation 
around the averages (triangles and error bars) are compared with experiment (filled circles and error bars). 


Global fit I Global fit II Global fit III 



1+2 3+4 5 1+2 3+4 5 1+2 3+4 5 


FIG. 2: Percentages of quark and glue components of /o(1370) obtained from Monte Carlo simnlation for global fits I (left), II 
(middle) and III (right) defined in subsection IV-A. In the first row, the components 1 to 5 on the horizontal axes respectively 
represent udud, {dsds + susu)/y/2, ss, (mm -I- dd)/\/2, and glneball G. In the second row, the percentages of total fonr-qnark 
(1-1-2), total quark-antiquark (3-1-4) and glue (5) are given. 
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Global fit I Global fit II 




Global fit III 



FIG. 3: Histograms of Monte Carlo simulations versus each component of /o(1370) obtained in global fits I (left set of 
figures), II (middle set of figures) and III (right set of figures). Each set consists of five figures, the vertical axes represent 
the percentages of Monte Carlo simulations, the horizontal axes (in each set) from left to right respectively represent the 
percentages of components udud, {dsds + susu) /%/2, ss, {uu + dd) /v^, and glueball G. Each component is divided into five 
percentage intervals (0-20%, 20%-40%, etc). For example, the first figure on the left shows that more than 90% of Monte Carlo 
simulations for global fit I, resulted in estimating the udud component of /o(1370) to be below 20%, and less than 10% of the 
simulations resulted in estimating the {dsds + susu)/'/2 component between 20%—40%, etc. 


where = 1.3. This leads to a subset S'ia C Si defined by 

Six = {p\p € Si : relation (IHTl) is upheld} . (33) 

The results for the components of /o(1370) over subset Six are shown in Fig. 0] (left panel) where the components at 
a given point in this subset are shown by a “-I-” together with the averages (triangles) and the standard deviations 
around the averages (error bars). We see that while zooming in on /o(1370) slightly shifts the components compared 
to the global fit I, the general structure remains the same, namely that the quark-antiquark components (particularly 
ss) remain dominant. 

Similarly, we search through set S’:! for a subset that better describes the overall properties of /o(1370) by imposing 
the condition 


Xii 


< (x 


II 1 


(34) 


with 


.,4370 _ ,4370 , ,4370 

Xii - Xm + Xr 


t/1370 

■X(r/r)> 


(35) 


where This leads to a subset Sux C S':: 


5'iia = {p\p € Sii : relation (IMll is upheld} . 


(36) 


The results for the components are shown in Fig. 0] (middle panel) which shows the same characteristics as those in 
Fig.m in which the quark-antiquark components (specially ss) dominate. 

Finally, we search through set for a subset in which the properties of /o(1370) are better described. We impose 
the condition 


xiY7° < (x 


III 1 


with 


,4370 _ ,4370 
XlII — Xm 


t/1370 

■ Xr > 


where = 0.9. This leads to the subset ^iha C S':!! 


'S'liiA = {p IP G <5111 : relation (l37t is upheld} . 


(37) 


(38) 


(39) 
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The results for the components are shown in Fig. |4] (right panel) which again shows a similar behavior as those in 
Fig. [5] where the quark-antiquark components (specially ss) dominate. 

This figure shows the robustness of the results and that when zooming in on /o(1370), while we get subsets that 
better describe the overall properties of this state, its components do not change much. As stated previously, this is 
not the case with other isosinglet states and in that sense /o(1370) is evidently a special case. Simulation histograms 
versus components are displayed in Fig. [5]showing the preeminence of quark-antiquark components in /o(1370). 


2. Zoom B 


We can further zoom in by applying a more stringent condition on each set and examine whether the components 
of /o(1370) retain their pattern observed above. For this purpose we search for a subset Sib G Si such that all decay 
ratios that were used as part of defining Si) are within their expected ranges (since the mass of /o(1370) has been 
fixed by WA102 collaboration to 1312 MeV, we do not impose that highly constraining codition in this zoom). This 
means that for points in Sib the decay ratios of /o(1370) should fall within their experimental values. 





■p3 -pS 

^ 7Z7T jKK 5 ^ r\r)/KK 


are within their experimental ranges 




(40) 


The results for the components are shown in Fig. [51 where again we see that the ss component is quite pronounced 
in agreement with the preceding discussions. We have given in Fig. [7] a comparison of the ss versus nfi components 
for the three global fits and the zooms discussed in this section. We see that for the majority of the simulations the 
ss component dominates. 

We impose similar strong constrains on the other global sets II and III but we do not find any subsets where all 
the inputs for /o(1370) can be met. 
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FIG. 4: Percentages of the quark and glue components of /o(1370) obtained obtained over snbsets Sia [relation (I33ll 1. Sua 
[relation (I36D [ and 5iii [relation (1391 . Components 1 to 5 respectively represent udud, [dsds + susu)/y/2, ss, (uu + dd)/y/2, 
and glue. In the second row, the percentages of total four-quark (1-1-2), total quark-antiquark (3-1-4) and glue (5) are given. 
The averages (triangles) and standard deviations (error bars) are shown. 
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Zoom lA Zoom IIA Zoom IMA 



FIG. 5: Histograms for Monte Carlo simulations versus each component of /o(1370) obtained over subsets Sia [relation (1331) : 
the left set of figures], Siia [relation (1361) : the middle set of figures] and 5iiia [relation (1391) : the right set of figures]. Each set 
consists of five figures, the vertical axes represent the percentages of Monte Carlo simulations, the horizontal axes (in each set) 
from left to right respectively represent the percentages of components udud, {dsds + susu) /%/2, ss, {uu + dd) / \/2, and glue. 
Each component is divided into five percentage intervals (0-20%, 20%-40%, etc). Eor example, the first figure on the left shows 
that more than 90% of Monte Carlo simulations in zoom lA resulted in estimating the udud component of /o(1370) to be below 
20 %, and less than 10% of the simulations resulted in estimating the {dsds + susu)/'/2 component between 20%—40%, etc. 


Zoom IB Zoom IB 



1 2 3 4 5 1+2 3+4 5 


FIG. 6: Percentages of the quark and glue components of /o(1370) obtained over subsets Sib [relation (I40|) ]. Components 1 to 
5 respectively represent udud, {dsds + susu)/ss, {uu + dd)/y/2, and glue. The percentages of total four-quark (l-|-2), total 
quark-antiquark (3-1-4) and glue (5) are also given (right figure). 


V. SUMMARY 


In this work we studied the internal substructure of /o(1370) within a framework that is designed to explore global 
properties of scalar mesons below and above 1 GeV based on various underlying mixings among two- and four-quarks 
as well as a glue component. The framework is based on a nonlinear chiral Lagrangian model with two scalar meson 
nonets (a four-quark nonet and a quark-antiquark nonet) and a scalar glueball. This framework has been used in series 


of prior works (refs. [74l| - 98j) on the subject and has given a coherent description of various low-energy experimental 
data. We investigated the 14d parameter space of the model by performing global fits to the properties of scalars. 
We determined sets of points that give an acceptable overall description of experimental data, and thereby computed 
the quark and glue components of /o(1370). It was observed that this state is dominantly a quark-antiquark state 
(particularly ss) with some remnant of four-quark and glue components. We tested the robustness of the conclusions 
made, by further zooming in on the properties of /o(1370) within the global context where the model has an overall 
agreement with experiment for all scalars below and above 1 GeV. It was observed that while individual components 
somewhat vary the overall features remains intact, namely that the ss component of /o(1370) is dominant. 

Since the Lagrangian for the 1 = 0 scalars (US studied in this work is constrained by the Lagrangian of the 
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25%1 


GFI GFII GFIII ZIA ZIIA ZIIIA ZIB 

FIG. 7: Comparing the ss and nh components of /o(1370). Vertical axis shows the percentage of simnlations in which the ss 
component is greater than the nh component. The horizontal axis gives from left to right the three global fits followed by their 
zooms A and zoom IB. 


I = 1/2,1 states studied in we test the stability of the results when the 7 = 1/2,1 Lagrangian parameters are 
relaxed and included in our global fit. Overall, we find that there is no noticeable change in the results and that 
the conclusion for the substructure of /o(1370) remains the same. The additional flexibility can be used to further 
investigate zoom B discussed in previous section. We find that global set III contains a point that describes all the 
experimental inputs for /o(1370) given in Table ITll i.e. 

S'liiB = {p\p & S'!!! : m[/o(1370)], are within their experimental ranges} . (41) 

For this point the components of /o(1370) are shown in Fig. |8]and further confirm the results found in this work for 
the substructure of this state. 

Finally, we point out the measurement of Bs —>■ J/^tt+tt” decay by the LHCb Il2,^ where it is reported that 
Bg —>■ J/'i/'/o(1370) is “firmly established.” Belle also has reported the same decay [l26j| . Since this decay proceeds 
through production of an ss pair shown in the schematic diagram of Fig. |9l we interpret this experimental result as 
some support for our prediction of a significant ss component in /o(1370) in the present study. 
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FIG. 8: Percentages of the quark and glne components of /o(1370) obtained over snbset Suib [relation (I41II ]. Components 1 to 
5 respectively represent udud, {dsds + susu)/y/2, ss, {uu + cid)/%/2, and glue. 
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/o(1370) 


FIG. 9: Production of /o(1370) in decay Bs —>■ J reported by LHCb [l25| . The decay proceeds via production of an 
ss pair and is consistent with the prediction of a large ss component in the substructure of /o(1370) presented in this work. 


Appendix A: Brief review of the mixing mechanism for 1 = 1/2 and 7=1 scalar states 

Although the scalar mesons above 1 GeV are expected to be close to the conventional quark-antiquark mesons [i| , 
a close look at their properties reveals that this expectation may not be completely supported by experiment. For 
example, if Aq(1430) and ao(1450) belong to the same quark-antiquark scalar meson nonet above 1 GeV, then it is 
rather surprising that 7^0(1430), an isodoublet with one strange quark is lighter than the isotriplet member ao(1450) 
that should not have any strange quarks in a pure quark-antiquark nonet. The masses of these two states are reported 
in in PDG Q: 


m [ao(1450)] = 1474 ± 19 MeV > m [AT/; (1430)] = 1425 ± 50 MeV (Al) 

which shows that the central value of the mass of isotriplet is about 50 MeV higher than the mass of isodublet. Even 
if we take the experimental uncertainties into account, which allows the masses to be comparable or get in the right 
order, still it does not make up for the quark model expectation in which the isodoublet is noticeably heavier than the 
isotriplet (for example, for the cases of the tensor and axial vector nonets, also p-wave nonets, the isodoublet is about 
100 MeV heavier than the isotriplet). Also some of the decay ratios of ao(1450) and Aq(1430) do not quite agree with 
the pattern that one would expect if these two states were members of a pure quark-antiquark nonet. These decay 
ratios from PDG [I| can be compared with the SU(3) predictions (given in parenthesis): P /P [Aq —>• ttK] = 

0.98±0.34 (1.51), P [ao ^ KK] /T [oq ^ nr]] = 0.88±0.23 (0.55) and P [oq ^ nr]'] /T [oq ^ nr]] = 0.35±0.16 (0.16). 
There are other similar deviations discussed in (^ . 

A natural question would be whether the deviations of experimental data for ao(1450) and 7^0(1430) from what 
is predicted if these two states were members of a pure quark-antiquark nonet can be understood based on a mixing 
of this nonet with the four-quark nonet below 1 GeV. This question was raised in and a mixing mechanism was 
put forward. This mechanism is based on a simple picture that provides a natural and effective description of the 
properties of 7 = 1/2,1 scalar mesons below and above 1 GeV (k(900), 7^q( 1430), ao(980) and ao(1450)) within in a 
nonlinear chiral Lagrangian framework (its extensions to the linear sigma model frameworks have been also studied 
in [s^-js^). The mechanism assumes that there are two scalar meson nonets around about 1 GeV (a lighter pure 
four-quark nonet N and a heavier pure quark-antiquark nonet N') and shows that allowing these two nonets to slightly 
mix with each other leads to a natural description of the properties of ao(1450) and 7irg(1430). The underlying reason 
that makes this mechanism successful is due to the reversed mass ordering of the “bare” (unmixed) states in the 
four-quark nonet N compared to those in the two-quark nonet N' shown in Fig. 1101 From the lightest to heaviest 
these four “bare” masses are as follows: The lightest “bare” mass is the 7 = 1/2 in nonet N which has one strange 
quark and is therefore lighter than the 7 = 1 of this four-quark nonet followed by the 7 = 1 state of nonet N' with 
no strange quarks and the heaviest 7=1/2 “bare” state in nonet N' which has one strange quark. Therefore, in the 
“bare” mass spectrum the two 7 = 1/2 states are farthest apart and the two 7 = 1 states are closest to each other. 
This reverse ordering turns out to be the magic behind this mechanism. To see how this works, we first remember the 
following simple property in small mixing e of two states of mass m 2 > mi that results in physical masses m 2 > mi 
described by mass matrices 


ml e 
e ml 


ml 0 
0 ml 


(A2) 
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We can easily show 


+ y (A3) 

where S = nn\ — m\. When the two “bare” masses are degenerate {mi=m2) the physical masses are m\ = m\ — e, 
+ e. This shows that when the two “bare” states of mass mi and m 2 mix, the physical masses split away 
from the “bare” masses and that the magnitude of this splitting is inversely proportional to the difference of the 
“bare” masses squared. Using this property we can see in Fig. [Ma) that when the two 1 = 1 states mix, since they 
are closer to each other, they split more than the two 7 = 1/2 states (since (5i < S 1 / 2 ), and as a result, a level crossing 
takes place where the ao(1450) get pushed above the 770(1430) and therefore we can understand the experimental 
data based on the mixing of a pure quark-antiquark nonet with a four-quark nonet. The case of = 0 is shown in 
Fig. linKb). Also shown in Fig. [TT] is an schematic quark line diagram for the mixing of two states with the same 
quantum numbers (one in the four-quark nonet N and another one in the two-quark nonet N') in which a possible 
rearrangement of quark lines to generate the two mixed physical states can be seen. In addition to describing the 
mass spectrum, the mixing also makes it possible to describe the decay ratios mentioned above. 

The Lagrangian density for the I = 1/2,1 scalars is developed in 

= -aTiiNN) - bTT{NNM) - a'Tr(iVW') - 6'Tr(WiV'7W) - yTr {NN '), (A4) 

+C'Tr(iV'a^</)Tr(a^</), (A5) 

where A4 = diag(l, 1, x) with x being the ratio of the strange to non-strange quark masses, and a, b, a', b', 7, A, C, 
A' and C are a priori unknown parameters fixed by experimental inputs on mass and decay properties of 7 = 1/2 
and 7 = 1 states below and above 1 GeV. 
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FIG. 10: Mixing mechanism for isodoublet and isotriplet scalar mesons below and above 1 GeV. Existence of a pure four-quark 
scalar meson nonet beneath a quark-antiquark nonet [Fig. (a)] results in a mass spectrum for 7 = 1/2, 1 states where the two 
isotriplets are closer to each other than the two isodoublets (i.e. < 51/2). This is due to the inverted mass spectrum in a 

four-quark nonet compared to that of a quark-antiquark nonet. Allowing the states with the same quantum numbers to mix 
leads to a level repulsion in which the isotriplets split more than the isodoublets and consequently a level crossing occurs that 
results in ao(1450) becoming heavier than the An (1430). This mechanism also works when the two isotriplets are degenerate 
in mass [Fig. (b)]. The mechanism developed in ref. |9^ . in the leading order, favored situation (b). The level repulsion also 
shows that the light scalar mesons below 1 GeV get pushed down and become lighter than expected. 
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Physical scalar 
above 1GeV 
(dominantly qq) 


Physical scalar 
below 1 GeV 
(dominantly 
9W9) 


FIG. 11: Schematic diagram for a possible rearrangement of the quark lines in the mixing process. 


Note that the last term in (IA4I) which induces the mixing between the two- and the four-quark nonets N' and N 
can be rewritten as 

- 7 Tr(iViV') = 

oc {QaQ’^) {qdq") = (eabc^V) (qd^) 

= (qdf) (qeq^) {qff) = ^ det [7V'(g)] . (A6) 

where N't(q ) = Qn. tf’- The determinant structure is similar to the contribution of instantons to the scalar channel 
discussed in (371 - 1^ . 

With 

"lao = 2(a + 5), 

= 2a + (1 + x)b, 

it is shown in [9^ that 

TOao = ma' = 1.24 GeV, to/Co = 1.06 GeV, = 1.31 GeV, 7 = 0.58 GeV^ (A8) 

The interaction coupling constants are also found [o^ from various decay widths of isodoublets and isotriplet states 
A = 1.19±0.16GeV\ A'=-3.37 ± 0.16 GeV\ (G = 1.05 ± 0.49 GeV\ C'=-6.87 ± 0.50 GeV\ (A9) 
Further details can be found in 


“K' 


= 2{a' + b'), 

= 2a' -b (1 + x)b', 


(A7) 
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